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Abstract: We study nonperturbative corrections to the transverse momentum dis- 
tribution of vector bosons in the Drell-Yan process. Factorizing out the Sudakov 
effects due to soft gluons we express their contribution to the distribution in the 
form of the vacuum averaged Wilson loop operator. We calculate the nonpertur- 
bative contribution to the Sudakov form factor using the expansion of the Wilson 
loop over vacuum fields supplemented with the expression for nonlocal gauge invari- 
ant field strength correlator. Although the Wilson loop is defined in an essentially 
Minkowski kinematics, the part of the nonperturbative contribution depending on 
the invariant mass of the produced vector bosons is governed by asymptotics of the 
correlator at large space-like (Euclidean) separations and therefore can be calcu- 
lated using conventional nonperturbative methods. Applying the results of lattice 
calculations we found that the obtained expression for the nonperturbative power 
corrections is in qualitative agreement with known phenomenological expressions at 
large transverse momenta and deviate from them at small transverse momenta. 
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1. Introduction 

Transverse momentum distributions of the vector bosons produced in the Drell-Yan 
process is one of the classical examples of successful QCD description of the phe- 
nomenology of hadron-hadron processes at high energies |l], ^, [|. Increasing accu- 
racy of the experimental data combined with the state-of-art resummed perturbative 
calculations allow to test perturbative QCD predictions and, at the same time, es- 
timate nonperturbative corrections to them associated with the confinement effects. 
It was proposed a long ago to parameterize the latter corrections by introducing ad- 
ditional phenomenological parameters into perturbative formulas 0, |3| . Their values 
were extracted from the comparison with the experimental data [|[]. Despite the 
fact that the obtained QCD predictions agree well with the experimental data || 
throughout a wide interval of invariant mass of the vector bosons, Q 2 , and their trans- 
verse momentum, q 2 , the understanding of nonperturbative effects in the Drell-Yan 
process is still missing. This should be compared with the situation with hadroniza- 
tion corrections to the deeply inelastic scattering (DIS) and inclusive observables in 
e + e~— annihilation. There, the operator product expansion (OPE) allows to organize 
the nonperturbative corrections in inverse powers of hard scale A2 n /Q 2n and identify 
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the corresponding nonperturbative scales A 2n as universal matrix element of compos- 
ite local operators in QCD like gluon condensate in e + e~— annihilation and higher 
twist operators in DIS. Since the OPE is not applicable to the Drell-Yan production, 
the physical interpretation of the measured nonperturbative scales remains unclear. 
Different approaches have been proposed to describe these scales within QCD || [/]]. 

In this paper we shall follow the approach developed in || and consider non- 
perturbative corrections the transverse momentum distribution d 2 a/dQ 2 dq 2 in the 
end-point region describing the production of the vector bosons with large invariant 
mass Q 2 and small transverse momentum q 2 <C Q 2 . It is well known that both 
perturbative and nonperturbative corrections are enhanced in this region. Indeed, 
perturbative QCD analysis shows that the cross section gets large perturbative cor- 
rections OL n s jq 2 log 2n_1 (Q 2 /g 2 ). These corrections are associated with multiple soft 
and collinear gluon emissions and they need to be resummed to all orders in the 
coupling constant |T[]. Applying the Collins-Soper-Sterman resummation formalism 
one gets || 

Q 2 dQ f d l dq 2 =or oE e 3 / (Sa e-^e- s ^^C jh ( Xl , Q, b)C jh ,(x 2 , Q, b)+Y(Q, q, x u x 2 ) . 

(1.1) 

Here, X\ = Q/\fse y and x 2 = Q/^/se~ y are momentum fractions of parton with 
relative rapidity y and the invariant energy of colliding hadrons s, (Tq = is 
the Born level partonic cross-section and Cjh are given by convolution of partonic 
distributions f a /h( x ) with perturbatively calculable coefficient functions Cj a 

C jh (x, Q, b) = £ f 1 ^c ja {x/t b, Q)f a /h(Z, b /b) (1.2) 

a Jx ? 

with 6 = 2e~ 7B and je being the Euler constant. Integration over two-dimensional 
impact parameter b ensures that the total transverse momentum of the vector boson 
equals q. In Eq. the large perturbative contributions due to soft gluons are 

factored out into the Sudakov form factor given by 

r Q 2 d h 2 ( (l 2 \ 

S(b 2 , Q 2 ) = j b2/b2 ^ [A(a s (k 2 )) log + B(a s (k 2 ))j (1.3) 

with the coefficient functions A(a s ) and B(a s ) known to two-loop order 0. The 
function Y(Q,q,Xi,x 2 ) describes the part of the cross-section (|1.1|) regular at q 2 <^ 
Q 2 . It is obtained by subtracting the logarithmically enhanced part from the fixed 
order result. 

The Sudakov form factor ( |1.3| ) resummes large perturbative logarithms a™ ln 2n ~ m (Q 2 b 2 ) 
and it suppresses the large b 2 , or equivalently small q 2 region of the distribution. Al- 
though it is defined for arbitrary b 2 it was realized that the perturbatively resummed 
formula ( |1.1| ) should be modified at large b 2 . 
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Firstly, the expression (|1.3|) was found by summing of all large logarithms log(b 2 Q 2 ) 
in the region Q 2 ^> 1/b 2 ^> Aq CD where perturbative QCD works. However, since in 
( |1 . 1|) the integration goes over all values of b 2 the modification of the resummation 
formula is needed for 1/b 2 ~ Aq CD . To this end it was proposed [[| to introduce 
the scale Qo below which the perturbation QCD fails to replace the impact parameter 
b by the function 

= b 
" " \l + b 2 QlY' 2 

in the expression for the Sudakov form factor, S, and for the distribution function, 
C, in ( p. . 1| ) . Then, at large b thus defined function b* approaches its maximum value 
b* = 1/Qo and perturbative expressions remain well defined. 

Secondly, we expect nonperturbative corrections become equally important at 
small transverse momentum q 2 , or large b 2 . The latter corrections appear suppressed 
by powers of both scales, Q 2 and q 2 . However, for small transverse momentum 
q 2 <C Q 2 it becomes reasonable to neglect all power corrections on the larger scale, 
Q 2 , and restrict consideration to nonperturbative effects on the smaller scale q 2 . The 
standard way to incorporate leading nonperturbative effects into ( |1 . 1| ) amounts to 
replacing the Sudakov form factor in ( |1 . 1| ) by the following expression [||, 0. 

S PT (b\ Q 2 ) - S PT ((6*) 2 , Q 2 ) + S nmiPT (b\ Q 2 ) (1-4) 
with nonperturbative part given by 

SnonPTO 2 ,*? 2 ) = 0(6) In tJ(- +<j>j/h(P,X 1 ) + (pj /h > (b, x 2 ) . (1.5) 

Here, and <f>j/h are some phenomenological functions whose form should be fitted 
to the experimental data. The simplest parameterization of the functions looks as 
follows 

(f)(b) = g 2 b 2 , <frj/h(b, xi) + (f)j/ h '{b, x 2 ) = gib 2 , (1.6) 

The values of the scales g\ and g 2 cannot be calculated in perturbative QCD and 
one finds them from comparison with the experimental data as ||, [| 

9l = 0.24 t° f 7 (GeV) 2 , g 2 = 0.34 t™8 (GeV) 2 . (1.7) 

There exist other proposals for the parameterization of nonperturbative effects in 
this process ||. 

In the present paper we study nonperturbative corrections to the resummation 
formula ( |1.1|) and calculate the power suppressing factor (|1.6|) using the notion of 
QCD vacuum condensates. The paper is organized as follows. In Sect. 2 we con- 
sider the origin of the Sudakov form factor in the resummation formula ( |1 . 1| ) . Using 
factorization properties of soft gluons we express this form factor as the expectation 
value of the Wilson loop operator built from soft gluon fields. In Sect. 3 we show 
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that the lowest order perturbative calculation of the Wilson loop expectation value 
reproduces known expression for the form factor Q1.3Q . At the same time, as shown 
in Sect. 4, to higher orders in perturbation theory this expectation value contains 
ambiguity associated with the contributions of infrared renormalons having the form 
of the power correction 0(b 2 A 2 ) ||. In Sect. 5 we calculate the nonperturbative 
power corrections to the Wilson loop by taking into account interaction of partons 
with vacuum gluon fields parameterized by nonlocal QCD condensates. In Sect. 6 we 
analyze the resummation formula for transverse momentum distribution with non- 
perturbative effects taking into account and compare it with the phenomenological 
ansatz (|1.5| )-( fL7D . Sect. 7 contains concluding remarks. 



2. Factorization of soft gluons 

The factorization formula ( |1.1| ) is valid up to corrections suppressed by a power of 
high scale Q 2 and it takes into account both perturbative and nonperturbative correc- 
tions on a smaller scale q 2 for q 2 -C Q 2 . The different factors entering the factorization 
formula ( |1 . 1| ) take into account contribution from different physical subprocesses - 
hard, collinear and soft - and have the following interpretation. All emissions of 
quarks and gluons with momenta collinear to incoming hadron momenta h and h! 
are factorized into two distribution functions f a /h(xi,b /b) and fb/h'(x2,bo/b). Hard 
emissions are described by the coefficient functions Cjh which, in general, are not uni- 
versal and depend on the produced vector boson. The soft subprocess describes both 
the initial state interaction between partons and final state radiation of soft gluons. 
It is factorized into the Sudakov form factor and will be the central object of our 
consideration. Its perturbative expression is given by ( |I~3l) and we anticipate that 
nonperturbative effects will modify it at the level of power corrections 0(Aqq D 6 2ti ) 
which we shall sum over all n. 

Let us demonstrate that the Sudakov form factor can be calculated in QCD as 
an expectation value of the Wilson loop operator built from soft gluon fields || . In- 
teraction of soft gluons with quarks can be treated using the eikonal approximation. 
In this way, one can obtain the perturbative expression for the Sudakov form factor 
by calculating the Feynman diagrams order by order in perturbation theory and re- 
placing quark propagators and quark-gluon vertices by their approximate eikonalized 
expressions. However, in order to go beyond the perturbation theory one should find 
an operator definition of the Sudakov form factor that does not refer to the Feynman 
diagrams and that is given entirely in terms of gluon and quark fields. Such definition 
exists and it is based on the Wilson line operators defined as 

W[C] = Pexp Ug J dZfAHzj) (2.1) 

with A^(z) being the gauge field describing soft gluons and C being an arbitrary 
integration path C in the Minkowski space. 



4 



The appearance of the Wilson line operators can be understood as follows |J. 
It is well known that in the eikonal approximation the interaction of quark with soft 
gluons leads to the appearance of the eikonal phase of the quark wave function 

-> %[-oo,x]^(x) (2.2) 



with $j,$p = 1. The phase oo,x] is given by the Wilson line operator (|2.1| ) 



evaluated along the classical trajectory of quark. For quark with momentum p^ this 
trajectory is a straight line that can be parameterize as = p^s + x M with s being 
a proper time 

$ p [— oo,x] = Pexp Ug J dsp fl A >M (ps + x)j . (2.3) 



Going over to the momentum representation, it is straightforward to verify that ( |2.2j ) 
and ( p.3| ) reproduce well known expression for the eikonalized amplitude of soft gluon 
emissions. 

In the case of the Drell-Yan process, the Sudakov form factor appears as a 
contribution of real and virtual soft gluons interacting with incoming quark and 
antiquark before their annihilation. Applying the eikonal approximation (|2.2|) we 
find that both quark and antiquark acquire eikonal phases which are factorized from 
the amplitude of the process into 

U BY (x) =T{$_ P2 [-oo,a;]$ Pl [:r,-oo]} , (2.4) 

where pi and p 1 ^ are momenta of quark and antiquark, respectively, x^ denotes the 
annihilation point and T stands for time ordering of gluon field operators. Projecting 
out the total eikonal phase (|2.4j ) onto all possible final states consisting of an arbitrary 
number of soft gluons, N > 0, with the total transverse momenta (—q) we can write 
the soft gluon contribution to the Drell-Yan cross-section as 

d 2 b 



a eik = Y,\(0\UM0)\N)\ 2 6^(k N + q)= J — e"^^ DY (6), (2.5) 
where 

WMb) = J2(^\UDY(0)\N)e- irkN (N\U j DY (0)\0) (2.6) 

N 

and \N) is the final state of iV soft gluon with the total momentum k 1 ^. Using 
the translation invariance property ?7dy(&) = e lPb [/DY(0)e _;iP6 with P^iV) = A;^| A^) 
being the total momentum operator, one can perform summation over all possible 
intermediate states to get || 

W DY (b) = (0|C/ DY (0)C/J Y (6)|0) = (0|Pexp [ig jT d^A"(z))|0) . (2.7) 

The integration contour Cdy entering this expression is shown in Fig. 1. It consists 
of two similar parts having the form of an (infinite) angle which are separated in 
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the transverse direction by the impact vector b. They belong to the plane defined 
by the quark momenta, pi and p%, and correspond to the eikonalized amplitude of 
the process and its conjugated counterpart. Two cusps located at the space-time 
points and b define the annihilation points. The cusp angle x is equal to the angle 
between quark momenta p\ and pi in Minkowski space 



coshx = {p\Pi)lp 2 



2p 2 



1 



(2.8) 



where Q 2 = (pi + p<i) 2 is the invariant mass of produced vector boson and we put 
p 2 to regularize light-cone singularity that appear for p 2 — > 0. We find 



Pi 



vl 



that for Q 2 /p 2 1 the cusp angle scales as x — hi and therefore the (kinematical) 
dependence of the Wilson loop on the cusp angle is translated into In Q 2 — dependence 
of the eikonalized cross-section (|2.5|) . 





(0,0,b) 



Figure 1: The integration path Cdy defining the eikonal phase (2/7). Straightforward 
lines correspond to classical trajectories of quark and antiquark with momenta pi and p2, 
respectively. 



We would like to notice that the Wilson loop in ( |2.7| ) should be averaged with 
respect to the state describing the spectator partons of both incoming hadrons. In 
this way one takes into account the possibility for quark and antiquark participating 
at the partonic subprocess to exchange by soft gluons with constituents of both 
hadrons before their annihilation. In what follows, we shall neglect such effects for 



the sake of simplicity and average (|2.7| ) with respect to the vacuum state. There 
are few reasons why such an approximation is reasonable. First, it is known [|§ 
that all perturbative effects of interaction between hadron constituents and particles 
participating in short distance partonic subprocess cancel in sum of all Feynman 
diagram up to corrections suppressed by a power of the hard scale Q 2 . Second, 
analysis of nonperturbative corrections to the transverse momentum distributions of 
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jets produced in e + e~— annihilation and deep inelastic scattering indicates [|IJ that 
the scale g 2 parameterizing nonperturbative corrections in Eqs. (|1.6|) and ( |1.7| ) is 
approximately universal and does not depend on the initial hadronic state. 

3. Perturbative contribution 

Using the expression (|2.7| ) for the Sudakov form factor we may try to take into 
account both perturbative and nonperturbative effects. The Sudakov form factor 
depends on Q 2 , lepton pair invariant mass of produced vector boson, and 6 2 , impact 
parameter conjugated to transverse momenta q. In the Wilson loop representation 
( |2T7|) , this dependence follows from the properties of the integration path Cdy in 
Fig. 1. In particular, the impact parameter 6 controls the transverse size of the 
integration contour. As a consequence, putting 6 = and using unitarity of the 
eikonal phase one finds that 

W DY (b = 0) = 1 (3.1) 

provided that all perturbative (light-cone and short-distance) singularities are ap- 
propriately regularized. This suggests to look for an expansion of the Wilson loop 
expectation value around 6 = 0. We expect that perturbative expansion will gener- 
ate In b 2 — corrections whereas nonperturbative effects will give rise to corrections in 
powers of b 2 . 

Let us start with perturbative calculation of (|2.7f) and show that the expression 
( |2.7|) reproduces the Sudakov form factor (2) in perturbative QCD at small b 2 . Using 
the definition ( [2.7D we get the following expression for the gauge invariant Wilson 
loop to the lowest order of perturbation theory: 

W (i) = I dx „ I dy v (Q\A»{x)A v {y)\Q) = l~Cc F [ dx» I dy u D^{x-y) 

(3.2) 

with Cf = (N 2 — 1)/(2N C ) and D^ u (x) being a free gluon propagator. To evaluate 
this integral we have to fix a gauge for a gluon field and, in addition, take care of 
possible divergencies of the integral. To this end we choose the Feynman gauge and 
introduce the dimensional regularization with D = 4 — 2e being the dimension of the 
x— space. Substituting the gluon propagator in the Feynman gauge and integrating 
it twice over the path Cdy we get the following expression for W^ 1 ' 

w» = 1 + sV -^/ I g 5 2^(* 2 )((^ - (^y) V E - 1), (3.3) 

where fi is a dimensionful parameter of the dimensional regularization, d D k = dk + dk-d D ~ 2 k 
with k± being light-cone components and k being (D — 2)— dimensional transverse 
momentum of gluon. Integration over k± leads to: 

p ik-b _ i 

{2n)D _ 2 p , (3.4) 
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where the cusp angle x was defined in ( |2.8|) . We notice that this integral does not 
have infrared divergences at small k 2 but it has ultraviolet divergence coming from 
large transverse momentum. Dimensional regularization provides a natural cut-off fi 
on maximal momenta of soft gluons and, as a consequence, the final expression for 
contains a pole in e 

w m = i _ ^±c F ( X coth x - i) r(1 ~ £) (7r/i 2 & 2 ) £ . (3.5) 

7r e 

Note that this divergence has nothing to do with conventional UV divergences in 
perturbative QCD since the residue of a pole in (|3.4| ) depends explicitly on the cusp 
angle \ which, by the definition ( j2.8| ), is a function of the kinematical invariants 
of the process. This divergence is introduced by the factorization procedure which 
separates the contribution of soft gluons from the exact amplitude of the process 
Ll] . Applying the eikonal approximation to the amplitude we correctly describe its 



IR behavior but change its UV asymptotics. Since by the definition only soft gluons 
contribute to the Wilson loop expectation value ( |2.7|) , we should put an upper bound 
on the value of the transverse momentum in ( [3.4|) . In the dimensional regularization 
this cut-off is given by fc 2 ^ = /i 2 . 

Specific UV divergences which we observed to one-loop order are general feature 
of Wilson lines. These divergences, the so-called cusp singularities, appear due to 
the fact that the integration path Cdy entering the definition Q2.7| ) of the Wilson 
lines contains cusps with angle \. The renormalization properties of the cusp diver- 
gences were studied in detail. It was shown that cusp singularities are renormalized 
multiplicatively and the corresponding anomalous dimension r cusp (o; s ,x) was found 
to the lowest orders in perturbative QCD JTT . 



Subtracting UV pole in the MS— scheme we get the one-loop renormalized Wil- 
son line as 

W 0) = 1 _ ^C F ( X coth X - 1) log^ (3.6) 

where the constant b was defined in (|1.3|) . Substituting expression ( |2.8j ) for the cusp 
angle x we encounter another singularity. Performing calculations of the partonic 
cross sections in perturbative QCD we neglect quark virtualities by putting p 2 = 
p\ = 0. In this limit the integration path in Fig. 1 goes along the light-cone and one 
expects additional singularities to appear. Indeed, using the definition fl2.8|) we find 
that the angle x is infinite at p 2 = and the one- loop expression for W^~> is divergent. 
As was shown in [11], this singularity originates from soft gluons propagating along 
the light-cone. The renormalization properties of light-like Wilson lines were studied 
in [12j . It was found that the renormalized Wilson loop having both light-cone and 



cusp singularities obeys the following renormalization group equation: 

d n . , d\dlnW(b 2 u 2 , 
% + (3{g) d- 9 d^Q 



^ + ^ gj ^^ --^M (3-7) 



8 



where r cusp (a s ) is the cusp anomalous dimension 

T cnsp (a s ) = ^C F + 0(a 2 ). (3.8) 

7T 

The general solution to this RG equation is given by || 

W(&V,tfV) = exp (-J^ ^-(r cusp (a s (k 2 ))log^ + T(a s (k 2 )) ] j^W (b 2 ,Q 2 ) 

(3.9) 

where T is the integration constant and Wo(b 2 ,Q 2 ) is the boundary value of the 
Wilson loop at the normalization scale /i 2 = b\jb 2 . Since /x has the meaning of 
the upper limit on the energy of soft gluons we can put jj 2 = Q 2 in ( |3.9| ). Then, 
the resulting expression for the Wilson loop expectation value ( |3"1^ ) is given by the 
product of two factors. One of them resummes all large perturbative logarithms in 
b 2 coming from integration over soft gluon transverse momenta above the IR scale 
1/b 2 . It coincides with the Sudakov form factor ( |1.3|) provided that the coefficient 
functions A(a s ) and B(a s ) are given by the anomalous dimensions, A(a s ) = T cusp (a s ) 
and B(a s ) = T(a s ). 



4. IR renormalons in Wilson loops 

The solution to the evolution equation ( pT9| ) is defined up to the factor W (b 2 ,Q 2 ) 
which gets contribution from soft gluons with transverse momentum below the scale 
1/b 2 and, therefore, cannot be calculated in perturbative QCD. It is this factor 
that provides nonperturbative contribution to the transverse momentum distribution, 
Eqs. ( |1.4| ) and ( P--7| ), and that we expect to be of the following form 

W (b 2 , Q 2 ) = exp (0(6) In ^- + O(ln Q)) . (4.1) 

Here, <fi(b) is some function of the impact parameter b 2 which can not be calculated 
perturbatively. In order to justify ( |4.1| ) and to determine the explicit form of the 
function <p(b) one has to perform a nonperturbative calculation of the Wilson loop 
expectation value ( |2.7| ). This will be done in the next Section. However, before going 
through an explicit calculations it becomes instructive to examine the general struc- 
ture of nonperturbative corrections to W (b 2 ,Q 2 ) by calculating the contribution of 
the IR renormalons to the Wilson loop expectation value (|2.7| ). The IR renormalons 
describe perturbative contribution to W \b 2 fi 2 , Q 2 / /i 2 ) coming from soft gluon with 
momenta k ~ Aqcd- The later turns out to be ambiguous at the level of power 
corrections ~ exp(— nj3o/a s (l/b 2 )) = (6 2 AQ CD ) n with the exponent n uniquely fixed 
by large order behaviour of the perturbative series. For the physical result ( |3.9| ) for 
the Wilson loop to be uniquely defined, this ambiguity should be absorbed into the 
nonperturbative power correction to the boundary condition Wo{b 2 ,Q 2 ). 
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To estimate the IR renormalon contribution to the Wilson loop we shall apply 
the renormalon technique and calculate (|2.7| ) in the so-called single renormalon chain 
(1— chain) approximation ||13|| . Applying the Borel transformation we can write the 
perturbative contribution to the Wilson loop as 



j^(l-chain) _ ^ 



du exp(— u/a s p )B[W] (u) 



(4.2) 



where the Born term was subtracted in the l.h.s. Calculating jy( 1 ~ cham ) we start with 
the lowest-order expression for the Wilson loop, replace the gluon propagator by the 
chain of fermionic bubbles and then apply naive nonabelianization to promote the 
quark contribution to the full one-loop expression for the (3— function. In this way, 



one finds the following representation for the Borel transform [13 



B[W}( 



u 



sin to 



l3o7rua s Jo 



A 



jy(l-cbain)p^ 



(4.3) 



where C is a renormalization scheme constant and C = —5/3 in the MS— scheme. 
Here, A 2 is the dispersive parameter which enters into perturbative expression for 
the Wilson loop in the form of a fictitious gluon mass. This allows us to write 
jy(i-chain)( A 2) in the form (pj) with $ + (k?) replaced by 5 + (k 2 - A 2 ). Integration 

over longitudinal components k± leads to the expression 



d V-2 k e tk-b _ 1 



(2tt; 



D-2 



k 2 + \ 2 



(4.4) 



which should be compared with (|3.4| ). Similar to (|3.6| ), it suffers from collinear 
singularities as p 2 — > 0, or x — * oo. To avoid these singularities it suffices to consider 
the following logarithmic derivative 



d 



d\nQ- 



-V^( 1_chain )(A 2 



(4.5) 



where K is the Bessel function of the second kind and D = 4 — 2s. Then, substituting 
731) into Q4.3D we obtain after some algebra 



B 



dW 
rflnQ 5 



(«) 



C f 



(iTb 2 ^ 2 y(b 2 fi 2 /4e c 



Til 



u — s) 



[u 



e)r(l+u)' 



(4.6) 



This expression for the Borel transform of the Wilson loop in the single chain approx- 
imation coincides with the similar expression obtained in || using different approach. 
Finally, substituting ( fOD into ( |4.2| ) and assuming that the contribution of a single 
renormalon chain exponentiates one arrives at 



dhxW 
dlnQ 2 



C h 



(7T/i b 



du T(l 



u 



U + S T(l + M 



(A 2 QCD 6 2 /4e c 



(4.7) 
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Let us analyse the Borel singularities of ( |4.6| ) and ( |4.7| ). We notice that at 
e — > the r.h.s. of ( [4.6|) contains a pole at u = 0. It has a UV origin and comes 
from integration over large transverse momentum in ( (4.4|) . As was explained in 



the previous section, the appearance of this pole is an artefact of the factorization 
procedure used to separate the contribution of soft gluon. Away from small u we 
put e = and observe that ( |4.6| ) and (|4.7|) acquires a series of IR renormalon poles 
located at integer positive u. This suggests to divide the integration region in (|4.7Q 
into two parts, < u < 1/2 and u > 1/2, and look for the solutions to (|4.7|) in the 
following form 

In W = In W {FT) + In W {lRR) (4.8) 

with pW PT ) getting a well-defined contribution from small u and W / ( IRR * ) taking into 
account IR renormalon poles. 

Integrating over small u in (|4.7| ) we expand the integrand in powers of u to get 



r/lniy( PT ) Ct? 00 (-Kii 2 b 2 e^ E Y r 1 / 2 



rfhiQ 2 "^9b^ C j ^ UttU y iv QCD 

Subtracting poles in MS— scheme and summing up the series we find 

d\nW^) _ C F Hble c /Al CD P) 

where we assumed that fi 2 , 1/b 2 ^> Aq CD . It is easy to check that the perturbative 
contribution to ( |3.9| ) satisfies this evolution equation and, as a consequence, 
coincides with the first factor in the r.h.s. of (13T 



Integrating in ( |4.7| ) for u > 1/2 we deform the integration contour around integer 
positive u = n and estimate the contribution of the IR renormalons to the Wilson 
loop as 

f/lnll/( IRR ) 00 ^ 1 

-^qt- = ]2/r^ QCD b 2 /biT - \r R (b) (4.io) 

with arbitrary coefficients c^ RR depending on the prescription that one uses integrat- 
ing the Borel poles and the thus defined IRR (6) being a function of b 2 regular at 
6 = 0. Comparing fl4.8|) with (|3.9| ) we notice that in order for the Wilson line expec- 



tation value W to be well-defined, the IR renormalon ambiguities of W^ 11111 ) should 
be absorbed into the definition of the boundary value Wo(b 2 , Q 2 ). The general form 
of the IR renormalon contribution can be obtained from ( |4.10| ) . It matches ( |4.1|) and 
IRR (6) provides IR renormalon contribution to the nonperturbative function <fi(b). 

5. Wilson loop in vacuum fields 

As was shown in the previous section, we expect nonperturbative effects to contribute 
to the power 0(b 2 A 2 ) corrections to the Sudakov form factor. To estimate these 
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corrections we shall use the representation of the Sudakov form factor as Wilson 
loop vacuum expectation value (|2.7|) . The calculation of Wilson loops in QCD has 
a long history. It is known that for small (Euclidean) size contours we may apply 
the perturbation theory to calculate the Wilson loop, while for large contours we 
expect essentially nonperturbative "area law" behavior. There were attempts to 



justify this behavior using different models of the QCD vacuum Trying to 

apply these approaches to the particular Wilson loop shown in Fig. 1 we notice the 
following important differences. First, the specific feature of the Wilson loop in Fig. 1 
is that it has essentially Minkowski light-cone geometry whereas the calculation of 
a Wilson loop is usually performed in the Euclidian QCD. Second, the integration 
contour in Fig. 1 is formally of an infinite size, but, at the same times, its definition 
involves a natural small parameter b 2 . At b = the integration path is shrunk into 
a point and the Wilson loop takes trivial value ( p.7| ). Then, for small b 2 <C 1/Aq CD 
the perturbation theory provides a meaningful approximation to W(b) whereas for 
b 2 ~ 1/Aq CD nonperturbative corrections become important. 

Following [|14| , |T5|| , we shall assume that for small b 2 nonperturbative corrections 
to W(b 2 ) come from interaction with vacuum fields. To evaluate the Wilson loop 
( |2~7l) we have to fix a gauge and then express the correlator of vacuum gauge fields 
in terms of gauge invariant quantities. This can be done most easily in the Fock- 
Schwinger gauge or fixed point gauge 

(x — x ) ■ A(x) = 0. 

The special feature of this gauge is that the gauge potential can be expressed in 
terms of field strength tensor as 

A^(x) = —(x - x ) u [ da aF^ v (x a + a(x - x )). (5.1) 
Jo 

Here, the fixed point xo is a free parameter of the gauge and the x$— dependence 
should cancel in gauge invariant quantities. To simplify the calculations we choose 

4 = ^ (5.2) 

with £ being an arbitrary gauge parameter. We will demonstrate that the final 
expression for the Wilson loop does not depend on £. 

Let us introduce the parameterization of the integration path of Fig. 1 as follows: 
CW = Ci + C n + C in + Civ-, where 

C I =p 1 s, C H = -p 2 t, C in = b-p 1 t, C IV = b + p 2 s, (5.3) 



as — oo < s < and < t < oo. Then, using the relation 05.1|)-(|5T3"D we evaluate the 
integral 



ig 



r rO rO rl 

/ dx-A(x)=ig / ds p^AJpi-s) = ig^ / ds daap^b v Fn V {pisa + b^{l — a)) 

JCi J-oo J-oo Jo 
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and changing the integration variables s — > —s/a and a — > 1 — a/£ we obtain the 
relation 

/ dx ■ A(x) = ig / / dapib u Fn U {—p\S + ba) (5.4) 
JGt jo jo 

in which £— dependence enters through the upper integration limit. Similar cal- 
culation of the integral along the part Cu , Cm and Civ leads to the following 
expressions 

fOO 

-ig 



ig 


/ dx 


A(x 




JCu 


ig 


1 dx 


A(x 




Cm 


ig 


/ dx 


A(x 




lc IV 



[ dt I dap^b" Fn U (— p 2 t + ba) 
Jo Jo 

poo pi 

J dt J dapiW F^ v (—pit + ba) 

poo pi 

J ds J dap 2 b u F fMU (—p2S + ba) (5.5) 

which depend separately on £. Combining together the expressions ( |5.4| ) and ( |5.5| ), 
we obtain the following relation 

p poo pi 

ig / dx ■ A(x) = ig / ds dab v (piFnJ—pis + ba) — p^F,^— p 2 s + ba)) (5.6) 
Jc Jo Jo 

in which the dependence of the gauge parameter £ disappears. 

Finally, substituting (|5.6| ) into (|3.2| ) we find that to the lowest order in a gauge 
field expansion the vacuum averaged Wilson loop is given by 



4iV 









<• 


g 2 






Jo Jo 







(5.7) 

Although this expression is £— independent, it is not gauge invariant. To restore 
the gauge invariance, two field strength tensors should be connected by a Wilson 
line to form a nonlocal correlator F® p (x)[P exp(ig Jy dz ■ A c (z)t c )] ab F^ x (y) with the 
group generators t a defined in the adjoint representation of the SU(N) group. Here, 
the integration path can be an arbitrary and we choose a straightforward line for 
simplicity. Any other choice leads to the expression for the Wilson loop which is 
different from ( |5.7|) by terms containing higher powers of gauge fields. Thus defined 
gauge invariant field strength correlator depends on the difference (x — y) and it can 
be decomposed onto two Lorentz tensors |i~6f : 

(0\gX P ^)[Pe i9 ^ dz - A{ %A(0)\0) = {g»„g pX - g, x g up )^(-x 2 + ie) (5.8) 

+ {g pa d^ - g m d p ){g\ a d v - sw<9 A )$ 2 (-x 2 + it) 

Here, $i and $2 are gauge invariant functions that receive both perturbative and 
nonperturbative corrections. In the special case of abelian gauge group the Bianchi 
identity e^^d^Fyp = implies that $i(— x 2 ) = 0. Although it is not the case for 
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QCD, this property leads to vanishing the lowest order perturbative QCD correction 
to $i. It is convenient to introduce the following representation 



roc , _ 

$i(-x 2 + it) = -i / dae~ iax ^ x (ia) 
J o 



d § 2 (-x + it 
d x 2 



: f°° dae- iax2 $ 2 {ia) , (5.9) 
J o 



where we took into account that eq. (|5.8| ) defines the function $ 2 up to an arbitrary 
constant and it is the derivative of this function that should enter into the final 
expressions for the Wilson loop. The correlator ( |5.8| ) receives both perturbative and 
nonperturbative corrections. To the lowest order of perturbation theory we use the 
definition ( |5.8| ) to calculate the functions $i(a) and $2(0") as 

$rv)=o, $r ) (-) = - i ^(iv 2 -i)/, 4 -^/ 2 - 1 . (5.10) 



Let us substitute the correlator (|5.8| ) into the expression ( |5.7| ) for the Wilson 
loop. Contracting the Lorentz indices, using the identities (pi ■ b) — (p 2 • b) — and 
performing integration over s in ( |5.7| ) one arrives at 



(P1P2) [ dsdte- iaipiS+P2t)2 = xcothx (5.11) 

Jo 2ia 

where \ is the angle between quark and antiquark momenta defined in (|2.8|) . Here, we 
put quark momenta off-shell, p\ = p 2 , 7^ 0, in order to avoid light-cone singularities. 
The remaining integral in ( |5.7| ) is given by 

fiiab 2 ) = tab 2 f 1 da ^e"^ 2 ^ 2 (5.12) 
Jo Jo 

and it can be expressed in terms of the error function erf (2) = Jq dt e - * 2 as 

r 1 2 
fix) = 2x da(l- a) e~ xa = e~ x - 1 + y/irx erlfy/x) (5.13) 
Jo 

for an arbitrary x. Finally, we obtain the following expression for the Wilson loop 

(E3 

W = l~ ^7(xcoth X - 1) / — \^{ia)f{iaP) - 2$ 2 (^)(1 - e^ 2 ) , (5.14) 
4iV Jo a z 1 J 



where the functions $1 and $2 parameterize the correlator ( |5.8|) . The following 
comments are in order. 

The dependence of the Wilson loop on the invariant mass Q 2 and impact param- 
eter b 2 is factorized in ( |5.14j ) into two different factors. The Q 2 — dependence is given 
by a "cusp factor" (%cothx — 1) which coincides with an analogous factor in the 
perturbative expression ( |3.4j ). At the same time, to evaluate the b 2 — dependence of 
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W we need to know nonperturbative functions $1 and $2- Notice that in the original 
expression for the Wilson loop (|5.7| ) the nonlocal gluon correlator was integrated over 
both space-time and like-like separations x 2 corresponding to the distance between 
different points on the contour Cdy, whereas the final expression (|5.14 ) contains the 
integral over parameter o and depends on the space-like vector b. This suggests to 
perform the Wick rotation: a — > —ia and rewrite ( |5.14j ) in the limit Q 2 — > oo as 

- / — <$> l {a)f{ab 2 )-2<$> 2 {<j){l- 



dlnQ 2 



AN 



00 da r ~ -, 

%(cr)f(ab 2 



a* 



(5.15) 



The advantage of this representation with respect to ( p.!4| ) is that the light-cone 
singularities cancel out in the r.h.s. of (5.15) and therefore it has a well defined limit 
as quark momenta go on-shell, pi = p\ = 0. 

The above mentioned properties of the Wilson loop (|5.14| ) and fl5.15| ) can be 
understood as follows. As was shown in the previous section, leading nonperturbative 
~ b 2 corrections to the Wilson loop come from the region of small transverse momenta 
whereas collinear logs ~ In Q 2 /p 2 appear from integration over small angles between 
gluon momentum and one of the quark momenta p% and p 2 - Since integration over 
small angles and small transverse momenta can be performed independently, the 
Q 2 — dependence of perturbative and nonperturbative expressions, Eqs. Q3.4D and 
( |5. 14| ) , respectively, coincide. The remaining integration over transverse momenta 
of gluon is translated into o— integration in ( |5. 14 ) with nonperturbative corrections 
absorbed into the functions $i(cr) and $2(0")- 

Expression ( |3.2|) defines first two terms in the expansion of the Wilson loop over 
vacuum gauge field. Natural question arises about the contribution of the remaining 
terms in this expansion. It can be estimated by applying the nonabelian exponenti- 



ation theorem |L7| which allows to distinguish higher order terms according to their 
color structure. It leads to the following general expression for the Wilson loop 



W = exp(u>i + W2 + 



(5.16) 



where w n is the total contribution to the Wilson loop W with "maximally non- 
abelian" color factors, w n ~ CpN™~ 1 , the so-called "webs". Expanding this expres- 
sion one can express the contribution to the Wilson loop with arbitrary color factor 
in terms of maximally nonabelian webs which can be defined without reference to a 
perturbative expansion. Comparing ( J5.16 ) with ( |5.14 ) we find that obtained expres- 
sion for W provides a lowest order contribution to the web W\. Then, applying the 
nonabelian exponentiation theorem we may exponentiate the lowest order expression 
(|5.14 ) and write the Wilson loop in the form ( |5.16| ) modulo corrections coming from 
higher webs. According to (|3.4| ), w\ has the following asymptotics in large Q 2 limit: 
Wi = 0(\aQ 2 ). It can be shown following ]TIJ, that this is a general property of 
the webs corresponding to the Wilson loop shown in Fig. 1, w n ~ InQ 2 . Therefore, 
being combined together in the exponent of ( |5.16 ) , higher order webs w n renormalize 
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the coefficient in front of collinear logarithm ~ InQ 2 in lowest web contribution w±. 
This property implies that the derivative d In W/d In Q 2 does not contain collinear 
logs and, therefore, it is Q 2 — independent as Q 2 — > oo in agreement with ( |5.15|) . 

Substituting ( |5.10| ) into (|5.15| ) and subtracting pole in the MS— scheme, one finds 
that perturbative contribution to the Wilson loop satisfies the following evolution 
equation 

d\nW^ T \Q 2 b 2 /b^_ _ rQ 2 dk 2 ± 

Jb*/b* 



cusp 



(5.17) 



dliaQ 2 hlitf k 2 

where r cusp given by ( |3.8| ). The higher order webs provide a higher order perturbative 



contribution to the anomalous dimension r cusp . Integrating (|5.17|) one reproduces 



perturbative contribution to ( |3.9| ). Applying the nonabelian exponentiation theorem 
to calculate nonperturbative contribution to the Wilson loop we arrive at the evo- 
lution equation ( |5.15| ) which is valid up to higher web contribution. The evolution 
equation (|5.15| ) allows to calculate the coefficient in front of In Q 2 but it does not fix 
the corresponding scale Ql 



W {b 2 ,Q 2 )=ex V 



AN 



log Ql 



da 



o a z 



$ 1 (cr)/(cr6 2 )-2<l> 2 (cT)(l-e 



(5.18) 

This expression was first obtained in [ 18 1 . It is valid up to higher order webs con- 
tribution and it is assumed that perturbative contribution is subtracted from the 
functions $i and $2- One finds similar expression for the Wilson loop using the 



stochastic model for the QCD vacuum [14]. There, it corresponds to the neglecting 
the contributions of higher cumulants of gauge strength fields. 



6. Asymptotic behavior of Wilson loop 

Let us consider the properties of nonperturbative contribution to the Wilson loop 
( 5.18 ). As follows from ( |5.18 ), the b 2 — asymptotics of the exponent is determined by 
the properties of the function f(ab 2 ) defined in ( |5.13| ). 

6.1 Small b 2 behavior 



Using the definition ( |5.9| ), we find that f(x) takes positive values for x > and for 
small x it has the following asymptotics 

2 3 4 

oc> lc oc 

f(x)=x 1 h..., for x < 1 . (6.1) 

JK J 6 30 168 ' V ; 

Thus, the exponent of the Wilson loop (|5.18|) admits the following expansion in 
powers of b 2 

Wo = exp (- In ^ [ Cl b 2 + c 2 b 4 + c 3 b 6 + 0(b s )] J , (6.2) 
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where the coefficients c n are defined as moments of the functions $i(cr) and $ 2 (cr): 

1 f°° da / ~ , . ~ , , 

c 2 



AN Jo (T 



24iV Jo 



j™ da($ x {(T) -6$ 2 (cx)) (6.3) 



1 roo 

c 3 = / da a $1(0") - 10<& 2 (<7 



The expansion (6^2) is in agreement with analysis of the IR renormalon ambiguities 
of perturbative series for W performed in Sect. 4. Dimensionful parameters c n are 
related to the nonlocal correlator (|5.8|) . The explicit form of these relations can be 
found from the definition (B.9I) as 



1 r°° r n 

c i = 7X7 / df2 N 5 *) + 2$ 2( f2 ) ( 6 - 4 ) 
4iV jo L J 

c 2 = [*i(0) + 6$' 2 (0)] (6.5) 

c 3 = - I ^[$ , i(0) + 10$' 2 / (0)], (6.6) 

where prime denotes differentiation with respect to x 2 . These equations imply that all 
dimensionful coefficients c n except of c\ are related to the short distance asymptotics 
of the nonlocal correlator (|5.8|) . Expanding the l.h.s. of ( |5.8|) in powers of x 2 one can 
calculate the first few terms of the small x 2 expansion of the functions $i(— x 2 ) and 
$ 2 (— x 2 ) in terms of local vacuum condensates as 

*i(-r) = ^P + ^^^ + 0(^ 4 ), (6.7) 

where (g 3 fF 3 ) = (g 3 f abc F° v F* p F^) , (g 4 j 2 ) = (g 4 j^) and j* is a quark current. 
Using these relations we evaluate the coefficients c 2 and C3 as 

c _ (9 2 F 2 ) _ 1 ( (9 3 fF 3 ) (g*j 2 }\ 

2 288N ' 3 1152iV V 5 3 / 1 ' ' 

We would like to stress that, in general, the scale C\ is not related to local vacuum 
condensates. 2 To find its value one needs to know the behavior of the correlator (|5l 
at large (Euclidean) distances. 



2 This is in agreement with the fact that there are no local gauge invariant operators in QCD of 
the dimension 2. 
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6.2 Large b 2 behavior 

The expansion ( |6.1| ) and ( |6.2| ) has a finite radius of convergence and it is not appli- 
cable at large x, or equivalently large impact vectors b 2 . Using the definition (|5. 
we find the asymptotic behaviour of the function f(x) at large x as 



f(x) = y/nx — 1 + e 
Its substitution into (|5.18[ ) yields 

w (b 2 ,Q 2 : 



1 

2x 



Ax 2 



+ ... 



for x ^> 1. 



cxp 



Q 



2 r 



<n Ql 



c \fb 2 + O(b 



where the dimension 1 coefficient Cq is defined as 



c 



(6.10) 



(6.11) 



(6.12) 



We notice that the function $ 2 (a; 2 ) does not contribute to Co and the large b 2 behavior 
of the Wilson loop is governed entirely by the function $i(x 2 ). 

Summarizing, we conclude that the Wilson loop has different behavior at small 
and large impact vectors b 2 : 



and 



lnW (b 2 ,Q 2 ) 



\nW (b 2 ,Q' 



In —2 CiO 

Qo 



— In Stt cqV 7 ^ 2 



Ql 



for 6 2 < A 2 



for b 2 > A 2 



(6.13) 



(6.14) 



where A is the characteristic length defined by the properties of the nonlocal correlator 
( |5.8|) . The nonperturbative scales cq and c\ are given by ( |6.12| ) and ( ]6.4|) , respectively. 
We recall that the expressions ( |5.18| ), (|6.13|) and (|6.14|) were obtained by solving the 
evolution equation ( 5.15|) and therefore they are valid up to corrections independent 
on Q 2 . In the expressions ( |1.5| ) and ( p..6| ) the later corrections are parameterized by 
the scale g\. 

6.3 Ansatz for the nonlocal correlator 

Calculating the scales c\ and cq we shall rely on particular nonperturbative ansatz 
for the nonlocal correlator ( |5.£| ). Since their values depend on the behavior of the 
nonlocal correlator (|5.8|) at large euclidian distances x 2 < 0, we may apply the results 
of lattice calculations of this correlator. 

is exponentially decreasing 



According to |T9| , [14| , the nonlocal correlator (J5T 
function of the space separation x 2 which can be parameterized at large distances as 



$i(f 2 



, _J£1 

$e * 



dx 2 



X 



$ e * 



(6.15) 
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with A being the correlation length, $ and $' being the overall normalization scales. 
This ansatz leads to the following expression for the function <&(cr) defined in (|5.8|) 



a ~ 3/2 exp 



4crA 5 



(6.16) 



and to similar expression for the function & 2 {°')- The lattice determination of these 
scales gives [[19] 

$ = 7.2 x 10 8 A|, A = 1/(183 Ai) (6.17) 
with the lattice QCD parameter 



A L = (0.005 ±0.0015) GeV. 



(6.18) 



In addition, |$'| « $/10 and therefore we may neglect the contribution of the 
$2 —function to ( |5.18| ). Then, substituting the ansatz (|6.15 ) into the definition ( |6.10| ) 
and ( p. 12 ) we get the following expressions for the coefficients 

a2 2$A 3 
Cl = — = 3584 A| , c = — ^ = 78 A L . (6.19) 



2V 



N 



Their values depend strongly on the lattice scale A^. Choosing this scale at the 
upper boundary A^ = 0.0065 GeV, or equivalently A = 0.17 fm, and using (16.171 ) we 
estimate the scales ( |6.19| ) as 



(latt) 



0.15 GeV 



(latt) 
c 



0.51 GeV. 



(6.20) 



Let us compare our predictions with the phenomenological parameterization Eqs. ( |1.5|) - 



1.7). To this end we rewrite ( |1.5| ) in equivalent form keeping the In Q— dependent 



part only 



lnWt Cn \b\Q 2 ) = -S nonPT = -\n^fb\ (6.21) 
We observe that this phenomenological expression matches . 13| ) at small b 2 as 



Q 2 92,2 



(phen) 



92 

2 



0.17 



+0.04 
-0.04 



(GeV) 5 



(6.22) 



This expression is in agreement with our calculation ( |6.20j ). 

It follows from our analysis that the parameterization ( |1.6| ) and ( |S.21| ) of nonper- 
turbative corrections to the Sudakov form factor is correct only for small b 2 A 2 , or 
equivalently for the transverse momenta q 2 ^> A -2 . We would like to stress that al- 
though the numerical expressions ( |6.17|) are sensitive to the particular ansatz for the 
nonlocal correlator ( |5.8|) , the general form of the small b 2 expansion (|6.2|) is uniquely 
fixed. According to ( |6.14| ), we expect that at large b 2 (or small transverse momenta 
q 2 <^ A^ 2 ) the asymptotic behaviour of the Wilson loop, and as a consequence non- 
perturbative correction to the Sudakov form factor, should have different form ( |6.14| ) 
with the scale Co given by ( |6.12| ) and ( |6.17| ). 
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2 



x=b/X 



6 



Figure 2: The profile function w(b/X) governing the b— dependence of nonperturbative 
contribution to the Sudakov form factor ( |6.23| ). The dash-dotted line describes the asymp- 
totic behaviour ( |6.14 ). 



To summarize, the nonperturbative corrections to small momentum distribution 
of lepton pairs in the Drell-Yan process are described by the expression ( 5.181 ) which 



is valid up to Q— independent corrections. The asymptotic behavior of the exponent 
of ( |5.18|) at large and small impact vectors is given by (|6.13|) and ( |6.14|) , respectively. 



Combining together ( |5.18 ) with nonperturbative ansatz for the nonlocal gauge field 



correlator, ( |6.15| ), one can find the expression for the nonperturbative factor that is 
valid for an arbitrary b 2 . Another advantage of the ansatz ( |6.16| ) that integration 
in ( 5.1 8| ) can be performed explicitly leading to the following expression for the 



$A_ 4 

Ql 



nonperturbative factor 

\nW (b 2 ,Q 2 ) = - In ^ -^r^( fo A), (6.23) 
with the profile function 

w (x) = e - x (3 + x) -3 + 2x (6.24) 

which is plotted in Fig. 2. Going over to the transverse momentum space we define 
nonperturbative weight function / n0 nPT(^ 2 , Q 2 ) as 

/„onP T (A: 2 , Q 2 )= j-^y e* S W (b 2 , Q 2 ) (6.25) 

Then, inserting the factor Wo(b 2 , Q 2 ) into the integrand of ( |1.1|) we find that nonper- 
turbative effects produce a smearing of the perturbative q— distribution, Eq. ( |1.1| ), 



20 



with the weight function f non pT(k 2 , Q 2 )- Substituting ( |5.18|) into ( |6.25| ) and taking 
into account ( |6.15|) we obtain the k— dependence of the weight function f n0 nPT{k 2 , Q 2 ) 
as shown in Fig. 3. 



0.08 




1 2 3 4 5 



Figure 3: The nonperturbative primordial distribution f n0 nPT(k 2 , Q 2 ) at Q = 91 GeV and 
Qo = 1.6 GeV. The dash-dotted line describes the Gaussian distribution corresponding to 

( TO - 

7. Conclusion 

In this paper we have studied nonperturbative corrections to the transverse mo- 
mentum distribution of vector bosons in the Drell-Yan process. Factorizing out the 
Sudakov effects due to soft gluons we were able to express their contribution to the 
distribution in the form of the vacuum averaged Wilson loop operator. We have 
demonstrated that the Wilson loop take into account both perturbative and non- 
perturbative corrections. The former are resummed through the evolution equation 
whereas the later define the boundary conditions for their solutions. Analyzing large 
order perturbative corrections to the Sudakov form factor we found that perturbation 
theory provides an ambiguous contribution to power corrections ~ (b 2 A 2 ) n associ- 
ated with the IR renormalons. The IR renormalon contribution suggests the general 
structure of nonperturbative corrections and it can be absorbed into nonperturbative 
boundary conditions for the Wilson loop. 

We calculated the nonperturbative contribution to the Sudakov form factor using 
the expansion of the Wilson loop over vacuum fields supplemented with the expres- 
sion for nonlocal gauge invariant field strength correlator. Although the Wilson loop 
is defined in an essentially Minkowski kinematics, the part of the nonperturbative 
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contribution depending on the invariant mass of the produced vector boson Q 2 is 
governed by asymptotics of the correlator at large space-like (Euclidean) separations 
and therefore can be calculated using conventional nonperturbative methods. Ap- 
plying the results of lattice calculations we found that the obtained expression for 
the nonperturbative power corrections is in qualitative agreement with known phe- 
nomenological expressions at large transverse momenta and deviate from them at 
small transverse momenta. 
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